
    GetRatio = function(theta, phi, eta){
    
        ratio = matrix(NA,5,2)
        rownames(ratio) = c("ak","lk","aK","lK","l0")
        colnames(ratio) = c("0","1")
        
        ratio["aK",] = theta
        ratio["lK",] = phi[1:2]
        
        denom.theta.l0 = (eta["10"]*phi["a1"]+1-eta["10"]) / 
            (eta["00"]*phi["a0"]+1-eta["00"])
        denom.theta.l1 = (eta["11"]*phi["a1"]+1-eta["11"]) / 
            (eta["01"]*phi["a0"]+1-eta["01"])
        denom.theta = c(denom.theta.l0, denom.theta.l1)
        ratio["ak",] = theta / denom.theta
        
        denom.phi   = (eta["01"]*phi["a0"]+1-eta["01"]) / 
            (eta["00"]*phi["a0"]+1-eta["00"])
        ratio["lk",] = phi[1:2] / denom.phi
        ratio["l0",1] = phi["a."] / denom.phi
        
        return(ratio)
        
    }